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Abstract 

We present an updated version of the FEWZ (Fully Exclusive W and Z production) 
code for the calculation of and 7* /Z production at next-to-next-to-leading order in 
the strong coupling. Several new features and observables are introduced, and an order- 
of-magnitude speed improvement over the performance of FEWZ 2.0 is demonstrated. 
New phenomenological results for production and comparisons with LHC data 
are presented, and used to illustrate the range of physics studies possible with the 
features of FEWZ 2.1. We demonstrate with an example the importance of directly 
comparing fiducial-region measurements with theoretical predictions, rather than first 
extrapolating them to the full phase space. 



1 Introduction 



Production of W bosons plays several important roles in hadron-collider physics studies. At 
the Tevatron, this channel furnishes the single most precise measurement of the VT-boson 
mass. The determination of the charge asymmetry between W + and W~ production at 
both the Tevatron and the Large Hadron Collider (LHC) provides important constraints on 
quark and anti-quark parton distribution functions (PDFs). High-energy production of W 
bosons acts as a background to searches for the new W bosons expected in many gauge 
extensions of the Standard Model, while production of the W in association with jets is a 
major background to searches for supersymmetry and other theories beyond the Standard 
Model with missing-energy signatures. Production of W bosons is copious at the LHC, and 
the statistical errors have already been rendered completely negligible with the few inverse 
femtobarns of data taken. Only systematic errors remain, arising from both experimental 
and theoretical sources. In many of the relevant observables, the experimental systematic 
errors are already at the few-percent level or smaller. 

The theoretical understanding of W production has reached an advanced stage. The 
inclusive 0{a 2 s ) QCD corrections to electroweak gauge boson production have been known 
for some time p]. Exclusive production, which is necessary for any realistic prediction or 
phenomenological study in a detector of finite acceptance, is technically challenging but has 
been achieved [21 El HI El El 13- The electroweak and QED corrections are known [U E], 
including the leading-logarithmic terms arising from multiple-photon emission [10J. Initial 
efforts toward the combination of higher-order QCD and electroweak corrections have been 
made [HI [T21 US], and even progress toward the exact mixed QED-QCD [TJ] and next-to- 
next-to-leading order QED corrections [15] has occurred. The remaining theoretical uncer- 
tainties are estimated to be at the percent level for the majority of relevant observables. 
These corrections are available to the experimental community in a variety of simulation 
codes. 

We have previously released a 2.0 version of the simulation code FEWZ (Fully Exclusive 
W and Z Production) that implements the NNLO QCD predictions for neutral-current 7* jZ 
production at hadron colliders, and allows for arbitrary kinematic cuts to be imposed [16J. 
FEWZ features a parallelized integration routine suited for running on modern computer 
clusters, allows for multiple, arbitrary kinematic variables to be binned during a single run, 
and automatically calculates PDF errors for the total cross section and all histogram bins. It 
has been used by collaborations at the Tevatron, the LHC, and RHIC in their physics stud- 
ies. The purpose of this work is to apply the advances in FEWZ 2.0 to the description of W 
production also, and to introduce further improvements to the FEWZ framework, resulting 
in a FEWZ 2.1 version. We have further improved the integration routine to provide an ad- 
ditional order-of-magnitude gain in speed. FEWZ 2.1 incorporates the LHAPDF format [T7] 
to allow all PDF sets of interest to be studied. Several new observables are introduced in 
the new version of the program. 

Our manuscript is organized as follows. We discuss the new features of FEWZ 2.1 in 
Section [21 Both the new available observables and the performance increase are detailed. In 
SectionElwe present phenomenological results for W physics at the LHC. We present compar- 
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isons of both integrated cross sections within the fiducial region and the charge asymmetry 
with ATLAS data. Several ratios of W + over W~ distributions are shown, both to illustrate 
the range of studies possible with the scripts distributed with FEWZ 2.1, and to highlight 
physics features of W production at the LHC The importance of directly comparing fiducial- 
region measurements to theoretical predictions, rather than first extrapolating them to the 
full phase space, is demonstrated using the W + over W~ cross section ratio. Only before 
the extrapolation is the measurement error small enough to allow for discrimination between 
different PDF sets. We conclude in Section HI 

2 New features 

Version 2.1 applies all improvements of FEWZ 2.0 p3] to W boson production. The user 
may now make either of two executables, f ewzw or f ewzz, which correspond to W or 7* jZ 
production, respectively. Most source code is shared between the two executables. The 
sector structure is slightly different owing to a few differing diagrams and symmetry factors 
at NNLO between W and 7*/Z production. Both executables take similar basic input files, 
and the same histogram input. The user of FEWZ 2.0 will not notice significant differences 
in compiling and running the new version of the program. We highlight below several new 
features of FEWZ 2.1. 

2.1 New observables 

We have added two new observables to further facilitate physics studies. 

• Beam thrust [H] has been added as an available histogram in FEWZ 2.1. Beam thrust 
is defined in FEWZ 2.1 as a sum over observed jets [19]: 

t b = ^Yl bVl eH%-Y| > (!) 

where Q and Y are the dilepton invariant mass and rapidity, and pk,T and % refer to 
the transverse momentum and pseudorapidity of the k-th jet. Beam thrust permits a 
central jet veto to be imposed without a jet algorithm, making it more easily amenable 
to resummation techniques. We include it to facilitate comparison between fixed-order 
results and resummation results, as well as for study in its own right. 

• The transverse mass has been added to allow for W production studies. Transverse 
mass is defined as 

M T = y 1 2p Tl E Trniss (l - cos{A<f)i tmiss )). (2) 

Pti refers to the transverse momentum of the charged lepton, Exmiss refers to the 
transverse momentum of the neutrino, and A(j)^ miss denotes the angular separation 
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between them in the transverse plane. Mj- has a Jacobian peak sensitive to the mass 
of the W without being dependent on the unobserved longitudinal momentum of the 
neutrino. It is defined analogously for the Z for comparison purposes. 

Both new observables are available for creating histograms; cuts on transverse mass are also 
accessible in the input file. 

2.2 New output features 

Several new possibilities for controlling the histogram output have been added. 

• The user now has the option to output cumulative histograms in addition to traditional 
ones. This can be useful to study the effect of cut placement on the variable of interest. 
This feature adds more than the convenience of avoiding adding histogram bins by 
hand: the resulting technical precision on a cumulative bin is better due to cancellations 
between bins. 

• Most experimental analyses use bin sizes that vary over the range of the histogram. 
Formerly this could be handled in FEWZ at the cost of multiple runs, but now the 
user has the option to input their own histogram bounds for any given histogram. 

• In addition to the scripts of FEWZ 2.0 which allow one to perform addition and 
multiplication operations on the output of separate runs, an asymmetry operator has 
been added which directly combines the results of W + and W~ runs in the following 

way: 

a w + — a w - f o\ 

A w = ■ [6) 

<Jw+ + a w- 

Both the statistical and PDF errors are propagated consistently. 

2.3 Speed improvements 

The sector combinations detailed in [16] have been reworked to provide a modest reduction 
in variance of the integrand, reducing the time required to reach the target precision. Some 
of these sectors have undergone "sector recomposition," which patches together the pieces 
resulting from sector decomposition [201 121] by undoing the variable transformations needed 
in the case of singularities in multiple endpoints of the phase-space integrals. This results in a 
more consistent mapping of the integration variables to parton phase space for similar sectors. 
There are now 133 sectors for f ewzz and 154 for f ewzw. Alternative techniques for reducing 
the number of required integrals when applying sector decomposition have been discussed in 
Ref. [22j. We have also been careful to eliminate redundancies caused by calling the same 
code with the same arguments repeatedly. Instead, the results of these calls (where they 
cannot be avoided) are stored in a cache or buffer. Finally, we have rearranged how certain 
variables are stored, so that access is more likely to occur contiguously in memory, resulting 
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in fewer CPU cache misses. This can make a large difference in the case of parton distribution 
functions, which are accessed frequently, or histograms, which are filled frequently. 

Altogether, the speed improvements are significant. In Fig.[TJ we reproduce Fig. 1 of [TB] , 
which shows integration precision as a function of time. FEWZ 2.1 has been added to the 
comparison. The cuts, histograms, PDF sets, and benchmark machine are the same as 
Ref. [16J, with the exception of the addition of the two new histograms described in Sec. 12.21 
The time required to reach a target precision is an order of magnitude less than needed for 
FEWZ 2.0. 

FEWZ Error vs. Time 
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Figure 1: Relative error versus time for different versions of FEWZ. 



3 Example results 

To demonstrate the capabilities of FEWZ 2.1, we present in this section several results 
relevant for the study of IV-boson production at the LHC We impose the following set of 
standard acceptance cuts on the W-decay products: 

M T > 40 GeV, E Tmiss > 25 GeV, 

PTl >20GeV, |7ft|<2.5. (4) 

We begin by studying the total cross section for W ± production subject to these cuts, 
together with the ratio of W + over W~ within this phase-space region. To illustrate the 
sensitivity of the cross section to different extractions of parton distribution functions, we 
show results for three NNLO PDF fits: those of HERA [23J, MSTW [21], and NNPDF [25] 
(to avoid clutter in the plots and tables we do not show the NNLO results of ABKM [26] or 
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JR [27J, which have been presented in our previous work on Z physics [16]). 



HERA: 


a w - 


= 2134i| pb; 


a w+ = 3265+H pb; 


a w + 


= 1.530 






a w - 




MSTW: 


a w - 


= 210611 Pb; 


a w+ = 3181+g pb; 


a w + 


= 1.510 






a w - 




NNPDF: 


a w - 


= 208511 Pb; 


a w+ = 3255l| pb; 


a w + 


= 1.563 








a w - 





+0.017 
-0.015 

+0.014 
-0.013 

+0.026 
-0.026- 



(5) 



The 68% CL. PDF errors are shown for each set. The parametric uncertainties arising 
from imprecise knowledge of the strong coupling and other parameters are not included in 
these results. Scale uncertainties have previously been shown to be subdominant to PDF 
uncertainties for electroweak gauge boson production at the LHC p2J [30] . We note that the 
technical errors arising from the integration error are less than 0.1%. The results for each 
PDF set are consistent, although the MSTW results for aw+ and the ratio are more than 
one sigma lower than the results for the other sets. We compare the results above with the 
ATLAS results for the aw+ / &w- ratio [28] : 

J = 1.542 ± 0.007(stat.) ± 0.012(sys.) ± 0.001(accept.). (6) 
a w- J fid 

We note that ATLAS combines both the electron and muon channels and extrapolates them 
to a common fiducial region when obtaining this number. The subscript 'fid.' denotes that 
the ratio is measured in this common region, which coincides exactly with the phase space 
defined by our standard acceptance cuts. The last error labeled accept, denotes the uncer- 
tainty contribution arising from extrapolation. The NNPDF and HERA PDF predictions 
are consistent within one sigma from the ATLAS measurement, while the MSTW result 
differs by over two sigma from the central value of the measurement. The situation becomes 
murkier if the ratio of inclusive cross sections is compared instead. In that case, the measured 
fiducial cross section is extrapolated to the full phase space using Monte Carlo simulation, 
introducing an additional source of uncertainty. The ATLAS measurement of the inclusive 
ratio is 



a w - 



1.454 ± 0.006(stat.) ± 0.012(sys.) ± 0.022(accept.). (7) 



Upon combining the various sources of uncertainty in quadrature, an error twice is large 
than the one in Eq. ([6]) is obtained. This result should be compared with the theoretical 
predictions for the inclusive cross section ratio: 



HERA 
NNPDF 
MSTW 



(%+/%-)inc = 1-4401S 
{<?w+/vw J ) inc = l-47llo;o2o; 
(a w+ /a w -) inc = 1.4351°;°^. 



All three sets now agree within the estimated one-sigma uncertainty with the results in 
Eq. ([7]). The full potential of the experimental measurement to distinguish between different 
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PDF extractions is only realized if the fiducial cross sections can be directly compared to 
theoretical predictions. With FEWZ 2.1, this becomes possible. We note that a comparison 
of CMS results [2H] with the W + over W~ inclusive cross section ratio obtained with different 
NNLO PDF sets was presented in Ref. [30], where a similar point regarding the importance 
of studying the fiducial cross sections was made. In that study, parametric uncertainties 
on the ratio of W + over W~ cross sections were also studied, and found to increase the 
PDF-only error by a negligible amount. We rely on this result to safely neglect them here. 

We next present results for several distributions, to illustrate both the functioning of 
FEWZ 2.1 and aspects of W physics at the LHC The transverse mass and rapidity distri- 
butions for W are shown in Fig. [2j The Jacobian peaks for Mt = My/ are clearly seen 
in both plots. Also visible is a similar, although less prominent, feature at Mt = 50 GeV 
caused by the cuts of Eq. (jlj). At leading order, the missing Et and the lepton momentum 
must be equal and back-to-back in the transverse plane, leading to the noted minimum. This 
constraint is relaxed at higher orders when the W ± can recoil against additional radiation, 
populating the region 40 GeV < Mt < 50 GeV. The stronger peaking of the up-quark PDF 
as compared to the down-quark distribution at x ~ 0.1 is visible in the enhancement of the 
W + rapidity distribution near \Y\ ~ 2. 

In order to illustrate new features present in FEWZ 2.1, we present in Fig. [3] results for 
the W~ beam thrust and the W + cumulative leading-jet pt distributions. No significant dis- 
tributions between the three PDF sets are found for either observable. Additional radiation 
in W production goes like CfC( s /tt. This is a smaller quantity than the Cao^s/tt relevant for 
gluon-initiated processes, and explains the peaking of the beam-thrust distribution at lower 
values than observed in Ref. [31] for Higgs production and the quick plateau observed in the 
leading-jet transverse momentum distribution. 

The scripts distributed with FEWZ allow ratios of kinematic distributions to be easily 
studied by combining the results of separate runs. We present below in Fig. H] several 
examples relevant to LHC phenomenology. The ratio W + over W~ production as a function 
of both lepton pseudorapidity and transverse momentum is shown in Fig. H] for the MSTW, 
NNPDF, and HERAPDF PDF sets. In the left panel, a separation between the sets is 
apparent at central rapidity, larger than the estimated PDF uncertainty. The ratio of lepton 
Pt distributions shown in the right panel reveals a harder distribution for the W~ than the 
W + . This can be understood by considering the leading-order kinematics. The pt and rj 
of the lepton can be written in terms of the center-of-momentum frame scattering angle as 
Pt — V^sin^ /2, rj = —In tan(#/2). The larger fraction of events at large rapidity for W + 
observed in Fig. [2] translates into a larger fraction of events at larger lepton pseudorapidity 
for W + than for W~ . The leading-order kinematics implies that W + events with large rj 
have scattering angles with small sin 6, and consequently smaller lepton pt than those for 
W~ production. 

The W charge asymmetry as a function of lepton pseudorapidity, defined as 

dA w _ da w + /dr/i — da w - /drii 
dr/i da w + / dru + da w - / dr\i 

is well-known to provide a strong constraint on PDFs. We compare in Fig. [5] the bin- 
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Figure 2: Bin-integrated cross sections for the W~ transverse mass (upper left panel), 
^"rapidity (upper right panel), W + transverse mass (lower left panel), and W /+ rapidity 
(lower right panel) for all three NNLO PDF sets. The standard acceptance cuts of Eq. 
have been implemented. The bands indicate the PDF uncertainties for each set. 



integrated results for Aw for the MSTW, NNPDF, and HERAPDF sets with recent results 
from ATLAS [28]. The HERAPDF and NNPDF predictions are in good agreement with the 
data over the experimentally available range of pseudorapidity, while the MSTW predictions 
are slightly lower throughout this range. 



4 Conclusions 

We have presented an update to the program FEWZ, version 2.1, which extends the func- 
tionality of FEWZ 2.0 to allow for the study of W physics at hadron colliders. Several new 
features have been added, including new observables, the possibility of cumulative histograms 
or ones with varying bin sizes, and an interface to the PDF library LHAPDF. Major speed 
improvements have been achieved; NNLO precision in the presence of standard acceptance 
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Figure 3: Bin-integrated cross sections for the W beam-thrust distribution (left panel) and 
the cumulative leading-jet px distribution in W + production (right panel). 
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Figure 4: Ratio of W 7 " 1 " over W production as a function of lepton pseudorapidity (left 
panel), and lepton p T (right panel). 



cuts can now be achieved on a single machine in hours. We have presented example runs of 
the new code for W observables with several NNLO PDF sets, both to show aspects of W 
physics at the LHC and to demonstrate the array of studies possible with the analysis scripts 
distributed with FEWZ 2.1. Comparisons to early LHC data indicate excellent agreement 
of data and theory as calculated by FEWZ. With FEWZ, the extrapolation of measured 
quantities to the full phase space and consequent increase in experimental uncertainty can 
be avoided. We look forward to the continued use of FEWZ in understanding the properties 
of electroweak gauge bosons at hadron colliders with unprecedented accuracy. 
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Figure 5: Charge asymmetry at NNLO for three different PDF sets, compared with recent 
data from ATLAS. 
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